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We summarize results from a study of spherically symmetric collapse of a 
Q^ ■ charged (complex) massless scalar- field ||]. We present an analytic argument 
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which conjecture the generalization of the mass-scaling relation and echoing 
phenomena, originally discovered by Choptuik, for the charged case. Fur- 
thermore, we study the behaviour of the self-similar critical solution under 
external perturbations - addition of a cosmological constant A and a charge- 
conjugation e. Finally, we study the scaling-relation of the black-hole charge. 
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\ Using an analytic argument we conjecture that black-holes of infinitesimal 
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Q ■ mass are neutral or obey the relation Qbh <C Mbh- We verify our predic- 
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Q\ . tions with numerical results. 
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I. INTRODUCTION 



^ . Gravitational collapse is one of the most interesting phenomena in general relativity. The 

dynamics of a spherically-symmetric massless scalar-field coupled to general relativity has 
two kinds of possible end states. Either the scalar field eventually dissipates away leaving 
spacetime flat or a black-hole forms. Numerical simulations of this model problem have 
revealed a very interesting phenomena - a kind of critical behaviour which is a feature 
of supercritical initial conditions very close to the critical case p = p* (p is a parameter 
which characterizes the strength of the initial configuration, and p* is the threshold value). 
More precisely, Choptuik found a power-law dependence of the black-hole mass on critical 
separation p — p* of the form: 

f p < p* 

I (p — p y p> p 
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Subsequently the same type of critical behaviour has been observed for other collapsing 
fields: the collapse of axisymmetric gravitational wave packets ||, the collapse of spherically 
symmetric radiative fluids [|J. In all these model problems the critical exponent (3 turned 
out to be close to the value originally found by Choptuik f3 ~ 0.37, suggesting a universal 
behaviour. However, Maison |J has shown that for fluid collapse models with an equation 
of state given by p = kp the critical exponent strongly depends on the parameter k. 

The second key feature of Choptuik's results is the universality of the precisely critical 
(p = p*) evolution at the threshold of black- hole formation - it was found that the critical 
solution has a discrete self-similar behaviour (discrete echoing with a period A). 

In this paper we summarize results from a study of spherically symmetric collapse of 
a charged (complex) massless scalar field. The main question we consider is whether it is 
possible to generalize Choptuik's results for the general charged case [Q. This generalization 
is not trivial since the introduction of charge- conjugation destroys the invariance of the 
evolution equations under the rescaling u — > au, r — > ar, where a is an arbitrary positive 
constant. This invariance characterizes the system in the special neutral case and is crucial 
for the self-similarity of the critical evolution. This lack of invariance hints that charge- 
conjugation might destroy the phenomena in the general situation. 

Using a semi-quantitative arguments, which is based upon different behaviour of the mass 
and charge under the rescaling r — > ar, we conjecture that the significance and influence of 
the charge decreases during the evolution and that the critical behaviour appears. Following 
this we provide a numerical evidence which confirms the generalization of the mass-scaling 
relation and echoing phenomena for the charged situation. 

So far, the mass-scaling relation of the black-hole, and the influence of perturbations on 
the critical evolution itself, have been studied in context of internal perturbations in the 
initial conditions, such as a deviation of the field's amplitude from the critical one. It is 
of interest to study the behaviour of the critical solution under external perturbations. We 
consider here the influence of the charge e as such a perturbation. We also consider the 
effect of another external parameter, the cosmological constant A. 

The plan of the paper is as follows. In Sec. II we describe the evolution equations. In 
Sec. Ill we describe the algorithm and numerical methods. In Sec. IV we describe our dis- 
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cretization and error analysis. In Sec. V we describe our theoretical predictions and compare 
them with our numerical results. Sections V.A and V.B establish both qualitatively and 
quantitatively as well, the generalization of the mass-scaling relation and echoing phenomena 
for the charged case. In Sec. V.C we study the behaviour of the critical evolution under two 
such external perturbations: addition of a cosmological-constant A to the critical solution 
and an addition of a charge-conjugation e to critical initial-conditions of a complex neutral 
scalar-field. We also study the possibility of forming black-holes from subcritical (p < p*) 
initial conditions using external perturbations. In Sec. V.D we study the charge-mass rela- 
tion for infinitesimal black-holes. We show that for p — p* — > the black-hole charge tends 
to zero more rapidly than its mass. From this we conclude that black-holes of infinitesimal 
mass, which can be created from near-critical evolutions, are neutral, or obey the relation 
Qbh <C Mbh- Our numerical results confirm this conjecture. We end in Sec. VI with a 
brief summary and conclusions. 



II. THE EVOLUTION EQUATIONS 

We consider a spherically symmetric charged scalar field <fi. This is a combination of 
two real scalar fields 0i,02, which are combined into a complex one = <f>% + ifo. The 
electromagnetic field is described by the potential A, which is defined up to the addition of 
a gradient of a scalar function. The electromagnetic field tensor F is defined as 2dA, i.e. 

F a b = 2A[6 io ]. 

The total Lagrangian of the scalar field and electromagnetic field is ||: 

L = ~((f>., a + ieA a <P)g a \^ b - ieA b <f>*) ~ -^-F ab F cd g ac g bd . (2) 

2 107T 

Where e is a constant and 0* is the complex conjugate of 0. 
Varying 0,0* and A a independently, one obtains 0: 

<P, ab g ab + ieA a g ab (2<P, b + ieA h <P) + ieA a , b g ab <P = , (3) 

and its complex conjugate, and: 

]-F ab ., c g bc - *e0(0* a - ieA a (jf) + 2e0*(0 ;a + ieA a 4>) = . (4) 
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We express the metric of a spherically symmetric spacetime in the form [0,|J: 



ds z = —g(u, r)g(u, r)du z — 2g(u, r)dudr + r 2 dVt 2 . 



(5) 



The radial coordinate r is a geometric quantity which directly measures proper surface area, 
and u is a retarded time null coordinate. Here dVt 2 is the two-sphere metric. 

Because of the spherical symmetry, only the radial electric field F 01 = —F 10 is nonvan- 
ishing. This choice satisfies Maxwell's equation: 



F\ah-A = . 



(6) 



We introduce the auxiliary field h: 



h = - hdr 



(7) 



Using the radial component of Eq. (f|), we express the charge contained within the sphere 
of radius r, at a retarded time u, as: 



and the potential as: 



Q(u, r) = Ante J r(h*h — hh*)dr . 
o 



A = J Qgdr 



+ —F ac F bd g cd + e 2 A a A b (f><j)* + Lg ab . 
An 



The nontrivial Einstein equations are: 



G , 



r g 



9>Tih )r h* r ; 



(8) 



(9) 



The energy-momentum tensor of the charged scalar- field is 0: 

Tab = -z (4>-a<P* b + 0* a 0;fe) + \ (-<f>;JeA b (f)* + <p* b ieA a <j) + <f)*JeA b <f) - <f). b ieA a <f)*^ 



(10) 



(11) 



G, 



l + r 9_ 9 

g g \g, 



(12) 
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Regularity at the origin requires g(u, 0) = g(u, 0). The boundary condition h(u, 0) = h(u, 0) 
forces us to integrate the equations outward, and impose the normalization g(u, 0) = 
g(u, 0) = 1, which corresponds to selecting the time coordinate as the proper time on the 
r = central world line. 

The solution at a given r depends only on the solution at r' < r. We integrate Eq. ( |TTD 
and obtain: 

" (h-h)(h* -h*) 



g{u,r) = exp 



Arc 



-dr 



(13) 



Using Eqs. ([□]) and (p~2|) , we obtain after integration: 

9(^ r ) = \ /f 1 "^) 9dr . 

o ^ ' 

In terms of the variable h, the wave-equation Eq. (|3]) takes the form 



(14) 



Dh = h u - ^gh r = ^- (g - g) (h - h) - ^ (h - h) g - % ^ r gh - iehA . ( I ~, ) 



2r 



Using the characteristic method, we convert the scalar-field evolution-equation Eq. ([15]) into 
a pair of coupled differential equations: 



dh 1 _ x /, tX Q 2 fl — \ icQ - 



2r 3 



2r 



gh — iehAo ; 



(16) 



dr 
du 



We solve these equations together with the integral equations (|7|), (|T3"|), (|I4|), 
The mass contained within the sphere of radius r at a retarded time u, is: 



M(u,r 

Using Eqs. (|H|) and (14), we express M as: 
M(u,r) = 



r_( I <T 
2 I g + r\ 



(17) 
and ©. 

(18) 



g 



dr + 



±91 

2 r 



(19) 
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III. ALGORITHM AND NUMERICAL METHODS. 



A numerical simulation of the special uncharged case was first performed by Goldwirth 
and Piran ||. Gundlach, Price and Pullin || used a version of this algorithm to study the 
scaling behavior of the mass of the black-hole for the special uncharged case. We have used 



a version of the algorithm of Refs. [p],|10| , for the neutral case, and we have generalized it 
for the charged case. However, the methods used in Refs |||9| are not accurate enough for a 
treatment of the critical solution itself, because each successive echo appears on spatial and 
temporal scales a factor e 2A ^ 31 finer than its predecessor. 

The main improvements of our version closely resembles that of Ref. [TIJ: Taylor expan- 
sion of the physical quantities near the origin, conservation of the number of grid points 
using interpolation and choosing the outermost grid point to be the ingoing light ray that 
hits the zero mass singularity of the critical solution (p = p*). 

To solve numerically Eqs. (|16D and ([17]) we define a radial grid r n , where n = 1, . . . , N. 
We should emphasize that h is a complex field: h = h\ + ih,2, and obtain a set of 3N coupled 
differential equations. 

= ^ (9n - 9n) (hm ~ hin) - |f (hm - hm) 9 



ri 



+ T^Qnhin + eh 2 nA 0n ; (20) 



n 



dh 



2/i 



1 - Q 2 

(9n — 9n) (h2n ~ ^2n) — T^IT (^2n — ^2n) 9n 

gjiin - eh ln A 0n ; (21) 



1 f-i \J It J It* / \ AIL - "All I r-^ 

du 2r n v J 2r: n 

2r, 

dr n 



du "2 s "' (22) 
where g and h satisfy the boundary conditions 

9i = 9i = 1 5 hu = /in ; h 21 = h 21 (23) 

The initial data for the Einstein-scalar-Maxwell equations is just the value of h on the 
initial data surface, u = 0, the value of the parameter e, and a numerical choice of the initial 
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position in r of each ingoing null lines of the grid. The algorithm proceeds as follows: first 
we integrate along r for a fixed u and find in succession the quantities h, g, Q, A and g. 

We use Eqs. (|7|), (|13D, (f|), @ and QUI) respectively in order to evaluate these quantities. 

The integration is carried out using a three-point Simpson method for unequally spaced 
abscissas (Even if the grid is evenly spaced initially, it will not remain so during the evolution 
||). We next use Eqs. (^0}|2^), to evolve h and r one time step forward. We solve the 3N 
ordinary differential equations using the fifth-order Runge-Kutta method [|11|] . This process 
is iterated as many times as necessary: i.e. until either the field disperses or a charged 
black-hole forms. 

The time step Au is determined so that in each step the change in r n is less half the 
distance between it and the null trajectory r„__i, i.e.: 

Au < T ' n ~ r "- 1 . (24) 

9n 

Once a null trajectory arrives at the origin r = 0, it bounces and disperses along u = const. 
to infinity. The grid-point is therefore lost when the light ray hits the origin. 
If for a given shell M > Q at some time then it is possible that 

^ (r , u) = 1 _!M^) + «!M (25 ) 

will vanish. We identify the formation of a black hole when there is an r value that satisfies 
t± = M(r±, u) ± (M 2 (r±, u) - Q 2 (r ±1 u)) 1/2 . (26) 

In this case r± are the horizons of the shell. 

As we approach the stage when a black-hole forms g n becomes infinite. In our algorithm 
this means that the step size Am — > 0. The numerical approach to the horizon is stopped 
eventually by an overflow of g or underflow of Au. We can, still, estimate where and when 
a black- hole horizon appears from the condition r + /r — > 1. 

In H| this method was used to study the scaling behavior of the mass of the uncharged 
black-hole. However, this method is not accurate enough for a treatment of the critical 
solution itself. The main improvements of our version for the charged case closely resembles 



that of Ref. [10] where this method was used to study the critical solution for the special 
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neutral case. We will now describe shortly the sources of inaccuracy and the methods that 
we use to overcome them. 

The first source of an inaccuracy arises from the fact that the expressions for the quan- 
tities h,g,g,Ao and M contains an explicit factor of 1/r which diverges to the origin. We 
overcome this by using a Taylor expansion of h in r: 

h = + h m r + /i (2 V 2 + 0(r 3 ) . (27) 

Then we expand the quantities h,Q, g,g, A and M in r, where the expansion coefficients 
are all functions of h,(°\ h^*, h^*, hP> and h^*. Thus, in order to treat the solution 
near the origin one needs to find only these coefficients. This is done by fitting the first 
three values of h to a second-order polynome: i.e. we solve equation for ri, r2 and to 
obtain h^°\h^ and h^. We then use these coefficients to evaluate the values of h,Q,g,g 
and Aq for the first two values of r. The values of these quantities for other values of r is 
then determined using a three-point Simpson method for unequally spaced abscissas. 

A second source of inaccuracy arises directly from the behavior of the critical solution 
itself - as the critical solution evolves each successive echo appears on spatial and temporal 
scales a factor e 2A ~ 31 finer than its predecessor. This problem is solved as follows: The 
number of grid points decreases during the evolution. Once a null trajectory arrives at the 
origin r = it bounces and disperses along u = const, to infinity. The grid point is therefore 
lost when the light ray hits the origin. As the number of points decrease by half we double 
the number of grid points and we interpolate the new grid points which are half way in 
between the old ones. We chose the outermost grid point in such a way that this ingoing 
light ray hits the zero-mass singularity of the critical solution itself. Thus the doubling of 
the grid gives us just the exact logarithmic scaling needed at the origin. 

IV. DISCRETIZATION AND ERROR ANALYSIS 

Our grid is highly nonuniform in u if a horizon forms - as we approach to the horizon 
the step size Am decreases rapidly. Furthermore, even if the grid is evenly spaced initially, it 
will not remain so during the evolution ||. However, when we consider two grids in which 
one has twice the number of grid points as the other, the coarser grid spacing will remain 
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twice the size of the finer grid, when we consider points at the same physical location. We 
will discuss, therefore, our numerical convergence in terms of the initial grid spacing. 

We denote the relative size of the initial grid spacing by I. From the condition ([24]) we see 
that the spacing Am is proportional to I. As we have mentioned, we treat Eqs. ( |16D and flTTp 
as ordinary differential equations in u, and we solve these using a fifth-order Runge-kutta 
method |TTJ, the error-term is 0(/ 6 ). 

The calculation of the quantities h, Q, g, g, A and M in these equations is nontrivial, and 
is given by Eqs. (|7|), (|S|), (|T3|), (|T^), ([]) and (TiTS| ) respectively. The integrals are discretized 
using three-point Simpson method, the error term in the integration is 0(l 5 f^), where 
is the fourth derivative of the integrand / ||TT] . 



As we go from one grid to a finer one (i.e. doubling the grid when the number of grid 
points reach half of the original number) we have to interpolate to obtain the values at the 
new grid points. The error term in the interpolation is 0(l A ). 

As was mentioned in |12||, the main source of error is the boundary conditions g = g 



and h = h at r = 0. We treat these boundary-conditions by approximating the true value 
of h(u,r = 0) using an interpolation according to Eq. (p7[), the error term is 0(l 3 ). The 
situation is even more complicated: The right hand-side of flip] ) contains an explicit factor 
of 1/r, which in the exact solution is canceled by the boundary conditions. 

Three reasons lead to the crucial importance of checking the behaviour of the solution 
near the origin, in order to establish our confidence in the numerical results: 

1. The risk of numerical instability caused by the explicit factor of 1/r. 

2. The solution at some value of r depends only on the solution at r' < r, so a numerical 
error in the quantity h near r = 0, would cause an error at each r, and in all the 
relevant physical quantities. 

3. The critical solution itself, which is the main issue of this work, appears on ever smaller 
spatial scales during its evolution, leading to the formation of a zero-mass singularity 
at r = 0, and the importance of our confidence in the numerical solution near the 
origin is therefore clear. 
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As we can see, h(u, r = 0) is the basic quantity which influence all the other physical 
quantities, and at each and every value of r. 

Figure 1 displays the error in h(u, r = 0). In this figure the initial data is of family (c), 
with amplitude A = 1.1, and e = 1 (see V) for a discussion of initial data). For these initial 
conditions the gravitational field is strong (A > A*), and the scalar-field undergo a terminal 
gravitational collapse into a charged black-hole. 
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FIG. 1. The convergence of h(u,r = 0) with decreasing grid size. The upper panel establishes 
visually the convergence of the code, by showing the real part of the scalar field hi (u, r = 0) for 
different relative grid spacing: 1/1600, 1/800, 1/400, 1/200, 1/100 and 1/50. The numerical 
errors are so small that the different six curves (for the different grids) actually overlap. The 
bottom panel establishes quantitatively the stability and convergence of the code. It displays the 
In of the error (compared to a reference solution with 1600 grid points) as a function of In of the 
initial grid size I. The slope is 3.03 ± 0.06. This slope indicates that the error near the origin is 
largely caused by our handling of the r = boundary conditions. The initial data is of family (c) 
with amplitude {A > A*) and e = 1. 

Due to the fact that there are no useful analytical solutions available, we use of the 
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(c) 

1/50-1/1600 



- (c) 

- Slope=3.035 



numerical solution itself as the reference solution. The reference solution was taken as the 
solution with 1600 grid-points, and it was compared with the solutions of 800, 400, 200, 100 
and 50 grid-points. 

The difference of h(u, r = 0) between a given calculation and the reference calculation 
was squared, summed (for different vales of u) and the square root was taken. The top part 
of Fig. 1 establishes visually that the code converges: this graph displays the real part of 
the scalar field hi(u,r = 0) for different relative spacing of the grid-points: 1/1600, 1/800, 
1/400, 1/200, 1/100 and 1/50. The numerical errors are so small that the six lines (for the 
different grids) actually overlap. 

The bottom part of Fig. 1 establishes empirically the stability and convergence of the 
code when we decrease the grid spacing: the numerical error varies as l m , where m = 
3.03 ± 0.06. The fact that the power is close to 3 at the origin is an indication that this 
error is caused by our handling of the r = boundary conditions. 

Furthermore, we have performed by a similar manner a calculation of the error in h(u, r = 
0.03). The top part of Fig. 2 displays again the fact that the numerical errors are so small 
that the different lines (for the different grids) actually overlap. The bottom part of Fig. 2 
establishes once again the stability and convergence of the code: the error varies as l m , where 
m = 4.01 ± 0.05. This value of the power indicates that in this regime the discretization 
error due to the interpolation which we use in order to maintain the number of grid-points 
is dominant. 
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FIG. 2. The convergence ofh(u,r = 0.03) with decreasing grid size I. The upper panel establishes 
visually the convergence of the code, by showing the real part of the scalar field h±(u, r = 0.03) for 
different relative grid spacing: 1/1600, 1/800, 1/400, 1/200, 1/100 and 1/50. The numerical 
errors are so small that the different six curves actually overlap. The bottom panel establishes 
quantitatively the stability and convergence of the code. It displays the In of the error ( compared to 
a reference solution with 1600 grid points) as a function of In of the initial grid size, I. The slope is 
4.01 ±0.05. This slope indicates that in this regime the discretization error due to the interpolation 
is dominant. The initial data is the same as for Fig. 1. 

V. THEORETICAL PREDICTIONS VS. NUMERICAL RESULTS 

In this section we present our theoretical predictions and numerical results for the grav- 
itational collapse of a scalar field, for both the uncharged and for the charged cases. The 
numerical results that we present arise from a study of several families of solutions, whose 
initial scalar field <f> = h profiles are listed in Table I: 
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Family 



Form of initial data 



(a) 

(b) 
(c) 
(d) 



Ar 2 exp 



/ V-0.2 V 
V 0.1 J 



/lexp(-44r) cos(lOOr) 



(r) 



Ar 2 exp 



r-0.25 
0.1 



+ iAr 2 exp 



( r-0.15 V 

V 0.1 ) 

Aexp(44r) cos(100r) + L4exp(-75r) cos(200r) 



Families (a) and (b) represent uncharged scalar-field while families (c) and (d) represent 
charged (complex) scalar-field. The amplitude A is the critical parameter p. 



A. THE CRITICAL SOLUTION 

In this section we discuss the critical solution (p — p*) itself, both for the previously 
studied uncharged case and for the newly studied charged one. First, we had to find the 
value of the critical parameter. This was done by a binary search until we found the value 
of the critical parameter to the desired accuracy. 

Let u* denote the value of u at which the singularity forms. We define 

T= - In [{u*-u)/u*\ ; (28) 



R = r/{u* -u) = {r/u*)e T . (29) 

In terms of these variables, the critical solution for the neutral case is characterized by 
discrete self-similarity ||, i.e. h(R,t) and other form- invariant quantity such as M/r or 
dMj dr are periodic functions of T. 

The generalization to the charged case is not trivial. Charge-conjugation destroys the 
invariance of the evolution-equations under the rescaling u au, r —> ar, where a is an 
arbitrary positive constant. This invariance characterizes the system in the neutral case and 
it is essential for the appearance of critical, self similar behaviour. The lack of this invariance 
raises the interesting question whether critical behaviour will appear in the charged case. 

Using a semi-qualitative argument, which is based upon different behaviour of the mass 
and charge under the rescaling r — > ar, we might be able to answer this question. As 
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the critical solution evolves its structure appears on ever smaller spatial (and temporal) 
scales. For the special neutral case we know that each successive echo appears on spatial 
and temporal scales which are a factor or 1 = e 2A « 31 smaller than its predecessor. Under 
the rescaling r — > ar we have: 

h^h; (30) 

9^9\ (31) 

Q^a 2 Q; (32) 

A -> aA ; (33) 

= <?o + (? e -> + a 2 g e ; (34) 

M = M + M e -> aM + a 3 M e , (35) 

where g , M are the parts of g and M which do not depend on e, and g e , M e are the additions 
in the charged case. Since a < 1 we learn from Eqs. (|30"1-|35|) that as the field evolves, the 
significance and influence of the charge on the evolution near the origin is reduced. 

Furthermore, as the oscilations proceed the scalar field approaches a real function times 



a constant phase [[13). This leads to an addtional decrease in the charge of the form: 

Q^a^Q; (36) 

A -> a^Ao ■ (37) 

9 = 9o + 9e -> 9o + a 2i 9e ! (38) 

M - M ■ M, — aM Q + a 2i M e , (39) 

where £ is a positive constant (see Section V.D). 

When both effects are combined together Q/M decreases with each echo, approximately 
from Q/M — > aSl + £)(Q/M). Looking at the right hand side of equation fll6|) we find that 
the various terms scale according to 

1 1 a (1+25) 

Yr [9 ~ {9o + 9e)} (h-h)^ —(g - g Q )(h - h) - ~^~Uh - h) ; (40) 

^h-k)g^ a —^-{h-h)g- (41) 
eQ a^eQ 

zF 9h -> ; (42) 

ehA -> a (1+5) e/M . (43) 
14 



From here we learn that as the evolution proceeds the last three terms in the right hand 
side of equation flTE| ) becomes smaller relative to the first term by a factor larger than a 2 
with each echo. From these arguments we expect to find that for near-critical evolutions the 
influence of the charge on the evolution should decrease with each echo. 

We conjecture that in the precisely critical case (p = p*) and in the limit of an infinite 
train of echoes, the influence of the charge on the evolution near origin is "washed out" and 
we expect to have the Choptuik solution. In practice, using the fact that a< 1 we expect 
the influence of the charge on the evolution to be negligible even after a small number of 
echoes. Thus, once echoing begins the solution will approach the neutral one rapidly. Our 
numerical solution verifies these predictions. Fig. 3 shows the quantity max(2M/r) as a 
function of T for near-critical evolution of families (a), (b), (c) and (d). The solutions of 
families (b)-(d) were shifted horizontally but not vertically with respect to family (a) in 
order that the first echo of each family will overlap the first echo of family (a). After an 
initial phase of evolution the quantity max(2M/r) settles down to a periodic behaviour in T. 
The period is A m 1.73 which corresponds to previous numerical results found by Choptuik 
for the neutral case. 



i , , | , i , | , , , r 

(a)+(b)+(c)+(d) 
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FIG. 3. Illustration of the conjectures universality of the critical evolution in the gravitational 
collapse of a charged (complex) scalar field. Max(2M/r) is plotted as a function of the logarithmic 
time T for the neutral families (a), (b) and for the charged families (c) and (d) (from near-critical 
evolutions) . The curves were shifted horizontally (but not vertically) in order to overlap the first 
echo of each family with the first one family (a). After an initial evolution the quantity Max(2M/r) 
settles down to a behaviour which is periodic in T. The period is A s=s 1.73, which corresponds to 
previous numerical results found by Choptuik for the neutral case. This Figure provides a numerical 
evidence for the conjectured generalization of the echoing behaviour for the charged case. 

Fig. 3 also provides a numerical evidence for the universality of the strong-field evolution 
of a critical configuration. In order to verify whether or not this universality exist for each 
and every value r we display in Fig. 4 profiles of Mj r as a function of R for u values at which 
the quantity max(2M/r) takes its maximum as a function of T. (It should be emphasized 
that Fig. 4 is composed of 16 profiles - 4 for each family.) The close similarity of the profiles 
illustrates two important properties: 

1. For each family by itself, the critical solution is periodic. Each successive echo appears 
on spatial and temporal scales a factor e 2A « 31 finer than its predecessor. 

2. The close similarity of the profiles illustrates the uniqueness of the critical solution, 
independent of the initial profile and charge. 
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FIG. 4. The profiles of M/r for each of the four families, are plotted as a function of the 
logarithmic coordinate R in those times at which the quantity Max{2M /r) reaches its maximum 
as a function of T . (It should be emphasized that this figure is composed from 16 profiles - 4 f or 
each family.) For each family by itself, each echo appears on spatial and temporal scales a factor 
e 2A w 31 finer than its predecessor. Furthermore, the close similarity of the profiles illustrates the 
uniqueness of the critical solution, independent of the critical profile and charge. 



We turn now to the power-law dependence of the black-hole mass. The critical solu- 
tion by itself does not yield the black-hole scaling relation; we should perturb the critical 
initial conditions, which has a self-similar character. This would lead to dynamical insta- 
bility - a growing deviation from the critical evolution towards either subcritical dissipation 
or supercritical charged black-hole formation. To describe the run-away from the critical 
evolution we consider a perturbation mode with a power-law dependence X(u* — u)~ a 
where A oc (p — p*). Assume the range of validity of the perturbation theory is restricted 
to some maximal deviation a from the critical evolution, i.e. the evolution is approximately 
self-similar until the deviation from the critical solution reach the value a. From here on the 
evolution is outside the scope of the perturbation theory - there is subcritical dissipation of 
the field or supercritical black-hole formation. In either case, the evolution from this stage 
on has no self-similar character. 

We assume the perturbation in the initial conditions develops into a charged black-hole 
the time u a (p — p*) requires in order to have a a deviation from the critical solution is given 
by the relation 



Of course, a bigger initial perturbation requires a smaller time u a - the horizon is formed 
sooner. The logarithmic time T a is given by 



B. SCALING BEHAVIOR OF THE BLACK-HOLE MASS. 




(44) 




(45) 



where depends on a and u*. 
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In a following paper |],|14| we prove that in order to have the total logarithmic time Tbh 
until the horizon formation one should add a periodic term F[ln(p — p*)] with a universal 
period, w, which depends on the previous universal parameters according to: w = A//3, so 
that the full dependence of Tbh on critical separation p — p*, for the neutral case as well as 
for the charged case, is given by 



Tbh = HP - P*) + FIMp - P*)] + h 



a 



(46) 



Fig. 5 displays —Tbh as a function of ln(a), where a = (A — A*) /A*, for the neutral 
families (a), (b) and for the charged families (c), (d) as well. The points are well fit by a 
straight line whose slope is 1/a ~ 0.37. This is consistent with the relation (f46|). 
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FIG. 5. The logarithmic time of a black-hole formation, plotted as —Tbh «s o, function of ln(a) 
(where a = {p — p*)/p* ) for the four families. The points are well fit by a straight line whose slope 



is 1/a f» 0.37. This is consistent with the relation (4-t). The oscillations seen above the straight 
line are physical [lj. 



The exponent a is related to (3 the critical exponent that describes the power-law de- 
pendence of the black-hole mass. We define as the mass after n echoes. Following this 
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definition we define 

M (o) = M (o) + M (o) =(i + C)M ( o) , (47) 

where M is the part of M which does not depend on e and M e is the charge contribution 
in the charged case. From and |39[) it follows that 

M (n) = M (o) e -nA + CM (o) e -(3+20nA _ ^g) 

We substitute (|45| ) into (f48|) and obtain 

\nM BH = c k + f3\n(p-p*) + 0[e 2 (p-p*) 2 ^ l+ V] , (49) 

and find that j3 = 1/a and c k is a family- dependent constant. Here we have assumed that 
M(°) and C do not depend on p — p*. 

Thus, according to the perturbation theory, the critical exponent (3 describes both the 
system's response to perturbations to the critical evolution, and the deviation rate from it. 
In a following paper [fl^JTf] []we show that one should add a periodic term \l/[ln(p — p*)} with 
a universal period w = A//3 to the scaling relation (pE9]). 

Fig. 6 displays ln(m) as a function of ln(a) for the neutral families (a), (b) and for the 
charged families (c) and (d) as well, where m is the normalized black-hole mass in units of the 
initial mass in the critical solution. The points are well fit by a straight line whose slope is 
(3 ~ 0.37. This is consistent with the prediction according to which (3 = 1/a. The measured 
value of the critical parameter f3 in the charged case agrees with the one of Choptuik for 
the special non-charged case. Thus, Fig. 6 presents a generalization of the mass-scaling 
phenomenon for the charged case. 



see also |L5|] 



19 



3 


i i i 1 i i i i 1 i i i i 1 
" (a) 


-2 


i i i 1 i i i i 1 i i i i 1 i i i i 1 

- m 




Slope=0.370 




- Slope=0.374 


4 












-4 




5 








6 




£ -6 




7 


, , I , , , , I , , , , I 


-8 





-20 -15 
ln(a) 



(o) 

Slops.0.371 



-25 -20 -15 -10 
ln(a) 



(d) 

Slopa=0.372 



-20 -15 -10 -5 
ln(a) 



-25 -20 -15 -10 
ln(a) 



FIG. 6. Scaling of the black-hole mass. ln(m) is plotted vs. ln(a) for the neutral families (a), 
(b) and for the charged families (c) and (d), where m is the normalized black-hole mass in units 
of the initial-mass in the critical solution. The points are well fit by a straight line whose slope 
is (3 ps 0.37. This is consistent with the predicted relation (5 = 1/a. The measured value of the 
critical exponent (3 in the charged case agrees with the one previously found for the special neutral 
case. Thus, this figure presents a generalization of the mass-scaling phenomena for the charged 
case. The oscillations above the straight line are physical [1,14,15]. 

In the neutral case the black-hole mass is proportional to its radius and either the mass 
or the radius could be the physical order-parameter. In the charged case the black hole 
radius equals M + (M 2 — Q 2 ) 1 ^ 2 and in general it is no longer proportional to M. However, 
as we approach the critical-solution, the significance of the charge decreases with each echo. 
The quantity Q/M becomes smaller as the black-hole mass gets smaller and it is impossible 
to distinguish between the two possibilities. 



C. EXTERNAL PERTURBATIONS TO THE CRITICAL SOLUTION 

So far, the mass-scaling relation of the black-hole, and the influence of perturbations 
on the critical evolution itself, have been studied in context of internal perturbations in 
the initial conditions, such as a deviation of the field's amplitude from the critical one. It 
is of interest to study the behaviour of the critical solution under external perturbations. 
In particular, we have studied the behaviour of the critical evolution under two external 
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perturbations: 



1. Addition of a cosmological-constant A to the critical solution itself. 

2. An addition of a charge-conjugation e to critical initial-conditions of a complex neutral 
scalar-field. 

Adding an external perturbation to the critical solution, which has a self-similar char- 
acter, is expected to yield dynamical instability - a growing deviation from the critical 
evolution toward either subcritical dissipation or supercritical black-hole formation. In the 
situation where the external perturbation to the critical solution (which, in the unperturbed 
case, is expected to form a zero-mass singularity) leads to a formation of a black-hole, we 
have studied the question whether or not there is a power-law dependence of the black-hole 
mass on the external parameter: A or e. 



The amplitude was first set to its critical value A* (for A = 0). We then add a 
cosmological-constant, A, and study its influence on the evolution of the critical solution. 
A > led to a dissipation of the critical solution and, on the other hand, A < was found 
to give a finite-mass black-hole formation. 

We have studied the behaviour of the black-hole mass as a function of the cosmological 
constant. This is analogue to the addition of an external magnetic field to a system of mag- 
netic moments and studying the magnetization dependence on the strength of the external 
field exactly at the critical temperature T c (at t — T c the magnetization is zero without 
an external magnetic field, just as the zero- mass singularity for A = A* and A = 0). The 
analogy arises from the fact that in both cases the external perturbation, magnetic-field or 
a cosmological-constant, forces the order-parameter to have a non-zero value at T = T c or 
at A = A*, correspondingly. 

For A ^ 0, the generalization of Eq. ( |HD is 



1. Addition of a cosmological constant A 



r 




(50) 



a 
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In addition, one should add the term 



A/2rg(h-h) 



(51) 



in the right-hand side of Eq. (|T^). Evidence for the universality and the power-law depen- 
dence of the black hole mass on the external parameter A is shown in Fig. 7 which displays 
ln(ra) as a function of In ( | A | ) for critical initial-conditions (A = A*). The points are well fit 
by a straight line, i.e. 



M BH oc |A| 



i/s 



(52) 



where 1/5 ~ (3. 
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FIG. 7. Power-law dependence of the black-hole mass on the cosmological constant. ln(m) is 
plotted vs. ln(|A|) for the four families. The points are well fit by a straight line whose slope is 1/5, 
where 1/5 ~ j3. This provides an evidence for the power-law dependence of the black-hole mass on 
external parameters. The initial conditions are for A = A* . 



2. Addition of a charge- conjugation e 

The amplitude was first set to its critical value A* for the families (c) and (d), with e = 0. 
For e = 0, the initial conditions represent uncharged complex scalar-field. We then add a 
charge-conjugation, e, and studied its influence on the evolution of the critical solution. We 
have studied the behaviour of the black-hole mass as a function of the charge-conjugation e 
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(keeping the amplitude on its critical value A*). The numerical results are shown in Fig. 8, 
which displays ln(m) as a function of ln(e 2 ). The points are well fit by a straight line, i.e. 

M BH oc |ep , (53) 

where 7 ~ 2(3. 
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FIG. 8. Power-law dependence of the black-hole mass on charge-conjugation e. ln(m) is plotted 
vs. ln(e 2 ) for the charged families (c) and (d). The points are well fit by a straight line whose slope 
is 7/2 ?3 (3. The initial conditions are for A = A* . 

In order to have a more detailed picture of the formation for a charged black-hole near the 
phase-transition we have studied the conditions required the formation of a charged black- 
hole and the dependence of its mass on the critical separation \p — p* \ and on the charge- 
conjugation e. This was done both for subcritical initial-conditions and for supercritical 
initial conditions as well, which are close to the phase-transition (\A — A*\ A*). 

Fig. 9 displays the dependence of the black-hole mass which forms from subcritical ini- 
tial conditions (A < A*) as a function of ln(a) and ln(e), for family (c). This provides a 
numerical evidence for the conjecture that external perturbation such as charge-conjugation 
can lead to a black-hole formation with subcritical initial-conditions (A < A*). Charge 
oppose gravitation but the electromagnetic energy density also contributes to the gravita- 
tional binding, and by doing so it permits the formation of a black-hole even from subcritical 
initial-conditions. 
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FIG. 9. The black-hole mass, ln(m), is plotted as a function o/ln(a) = ]n(A* — A)/A* andln(e), 
for subcritical initial conditions (A < A*). The charge opposes gravitation but the electromagnetic 
energy also contributes to the gravitational binding, and by doing so it permits the formation of a 
black-hole even from subcritical initial conditions. Black holes do not form in the flat regime. A, 
B, C and D are normalization constants: A = 1, B = 16, C = -1, D =-11.061 

The critical charge-conjugation e*, needed to obtain large enough deviations from the 
otherwise evolution, followed by black-hole formation from subcritical initial conditions, is 
well described by a power-law dependence on critical-separation \A — A*\: 



e* oc (p* — pY 



(54) 



where e = 0.5 (see Fig. 10). 



FIG. 10. Power-law dependence of the critical charge- conjugation e* on the critical separation 
\A — A*\, for subcritical (A < A*) initial conditions, e* is the minimal charge- conjugation needed 
for black-hole formation from subcritical initial- conditions. 
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Fig. 11 displays the dependence of the black- hole mass formed from supercritical initial 
conditions (A > A*) as a function of ln(a) and ln(e), for family (c). We learn that near 
the phase-transition, the larger the critical separation p — p* and the charge- conjugation e, 
the larger is the black-hole mass. The slope of the right edge of the surface is (3, which is 
consistent with Eq. The slope of the left edge of the surface is 2/3, which is consistent 

with Eq. (pi). 




FIG. 11. The black-hole mass, ln(m), is plotted as a function of ln(a) = ln(^4 — A*)/ A* and 
ln(e) ; for supercritical initial conditions (A > ^4*). Near the phase-transition, the larger the critical 
separation p — p* and the charge- conjugation e, the larger is the black-hole mass. The slope of the 



right edge of the surface is (3, which is consistent with Eq. (4&)- The slope of the left edge of the 
surface is 2/3, which is consistent with Eq. (55). A, B, C, D are normalization constants: A = 1, 
B = 8, C = 1, D = 23.061. 

It should be emphasized that the tight relation between the various critical-exponents 
/3, 5 and 7, which all describe the instability of the critical evolution under a variety of 
different perturbations, is a strong evidence supporting the conjecture that there exists one 
mechanism which can explain the power-law dependence of the black-hole mass on the 
various parameters, both for internal perturbations in the initial conditions, and for external 
perturbations as well. 
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D. THE CHARGE-MASS RELATION FOR BLACK HOLES 



We have shown in 5. A that the black-hole charge tends to zero faster than its mass for 
p — p* — > 0. We demonstrate now that Q — > as a power law with a critical exponent 
larger than 2(3 (see also [13] for an independent analysis of this problem).. Define as 



the charge after n echoes. From and |36]) it follows that 

QW = Q(°) e -( 2 +?)« A . (55) 
Substituting ( [45|) into (|55|) and assuming that Q'- -' does not depend on p — p* we obtain 

ln|Q B ff| = (2 + 0/3HP~P*) + dk , (56) 

where dk is a family-dependent constant. 

From this analytic argument we deduce two conclusions: 

1. The black- hole charge is expected to have a power- law dependence on critical separa- 
tion p — p*, where the critical-exponent rj is closely related to the critical-exponent /3 
of the mass according to r\ — (2 + 

2. The black- hole charge tends to zero with p — p* more rapidly than its mass. 

Fig. 12 displays lng^ as a function of ln(a) for near-critical black-holes, where qsH is 
the normalized black-hole charge in units of the initial- charge in the critical-solution. The 
points are well fit by a straight line whose slope is rj ~ 0.88. This is consistent with the 
relation (|56T). 
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FIG. 12. Illustration of the conjectured charge- scaling relation\5^. In(qBii) is plotted vs. ln(a) for 
near-critical black-holes, where qsH is the normalized black-hole charge in units of the initial- charge 
in the critical solution. Data from the charged family (c) is shown. The points are well fit by a 
straight line whose slope rj obeys the relation rj > 2(3. Thus for p — p* — > the black-hole charge 
tends to zero more rapidly than its mass. 

Another way to determine the value of r\ is directly from Eqs. (|55] and ^). Fig. 13 
displays lnq( n > as a function of n, the number of echoes, along the critical solution. The 
slope is — (2 + £) A w —4.133. Using A ps 1.73 and (3 ~ 0.37 we find 77 ~ 0.883, in agreement 



with the prediction, 77 = 0.883 ± 0.007, of Gundlach and Martin- Garcia [13 . 
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FIG. 13. The decrease of the charge with each echo during the critical evolution. The points are 
will fit by a straight line whose slope is — (2 + £)A s=s —4.133. 

The data shown in Figs. 12 and 13 come from the charged family (c). For family (d) 
we have found that in general the charge increases as the critical separation p — p* increases 
although this correlation is not well described by a power-law dependence. The reason for 
this situation probably comes from the initial stage of the evolution (note that the initial 
data contains layers of both positive and negative charge whose relative magnitude depends 
on p. Consequently, the assumption that does not depend on p — p* breaks down for 
this configuration. 

Despite of this, our numerical results qualitatively confirm our prediction - the black-hole 
charge tends to zero with p—p* faster than its mass. From here we conclude black- holes with 
infinitesimal mass, which, according to this mechanism, can be created from near-critical 
evolutions, are neutral, or obey the relation Qbh <^ M BH . 
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VI. SUMMARY AND CONCLUSION 



We have studied the spherical gravitational collapse of a charged (complex) scalar-field. 
The main issue considered is the generalization of the critical behaviour, originally discovered 
by Choptuik for neutral fields, for the general charged case. This generalization is not 
trivial, since the introduction of a charge-conjugation destroys the invariance of the evolution 
equations under the rescaling u — > au, r — > ar, which characterizes the system in the neutral 
case and is crucial for the self-similarity of the critical evolution. This lack of invariance 
hints that charge-conjugation might destroy the phenomena in the general situation. 

However, we have shown that the significance and influence of the charge decreases during 
the evolution and consequently the critical behaviour appears. As p — p* — > the black-hole 
charge tends to zero faster than its mass, i.e. Qbh/Mbh — ► 0. Thus, we conjecture that 
black-holes of infinitesimal mass, which can be created from near-critical evolutions, are 
neutral, or obey the relation Qbh <C M BH . Consequently, we find both the mass scaling 
relations for supercritical solutions and the echoing phenomenon for the critical solution. 
The charge of the black hole also depends, as a power law on the separation from criticality. 
The critical exponent, is rj > 2f3. 

We have also studied the response of the critical solution to external perturbations. In 
particular, we have studied the behaviour of the critical evolution under an addition of a 
cosmological-constant A, and an addition of charge-conjugation e to critical initial-conditions 
of a complex neutral scalar-field. We have found a power-law dependence of the black-hole 
mass on external parameters A and e. 

Under the influence of these external perturbations it is possible to form black-holes 
from subcritical (p < p*) initial conditions. Once more, the critical charge-conjugation e*, 
which is needed to form black-holes from subcritical initial-conditions, is well described by 
a power-law dependence on critical separation p* — p. 
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